Charge and spin criticality for the continuous Mott transition in a two-dimensional 

organic conductor 
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We study the continuous bandwidth-controlled Mott transition in the two-dimensional single- 
band Hubbard model with a focus on the critical scaling behavior of charge and spin degrees of 
freedom. Using plaquette cluster dynamical mean-field theory, we find charge and spin criticality 
consistent with experimental results for organic conductors. In particular, the charge degree of 
freedom measured via the local density of states at the Fermi level shows a smoother transition than 
expected for the Ising universality class and in single-site dynamical mean-field theory, revealing 
the importance of short-ranged nonlocal correlations in two spatial dimensions. The spin criticality 
measured via the local spin susceptibility agrees quantitatively with nuclear magnetic resonance 
measurements of the spin-lattice relaxation rate. 
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The Mott metal-insulator transition (MIT) is a 
paradigmatic example of a correlation-induced phase 
transitionii Its physics is generically contained in the 
single-band Hubbard model, which is parametrized by 
the local Coulomb repulsion U, the bare bandwidth W, 
and the average electron density n. Two MITs are dis- 
tinguished: first, the bandwidth-controlled Mott tran- 
sition at fixed filling, where an insulator turns into a 
metal by increasing W/U, typically realized in experi- 
ments through chemical or hydrostatic pressure, and sec- 
ond the filling-controlled Mott transition at fixed U/W, 
where the system becomes metallic upon adding elec- 
trons or holes. A paramagnetic Mott transition is often 
superceded by antiferromagnetic ordering unless the sys- 
tem is frustrated or the temperature high enough so that 
the bandwidth-controlled MIT proceeds from a param- 
agnetic metal to a paramagnetic insulator. Typically, the 
paramagnetic MIT is a discontinuous first-order transi- 
tion at low temperatures with a continuous critical end 
point. 

The scaling behavior at the critical end point, or in 
short Mott criticality, for the bandwidth-controlled MIT 
has been probed experimentally for the charge degree of 
freedom via the dc conductivity, e.g. for (Vi-xCtx)20^ 
or the quasi-two-dimensional k-(ET)2Cu[N(CN)2]C1 (ab- 
breviated as «;-Cl)j^>^ The latter belongs to a class of 
layered organic charge transfer salts,— >^ which are both 
low-dimensional and geometrically frustrated, such that 
magnetic order is suppressed at the temperature of the 
critical end point of the Mott transition. These organic 
salts therefore allow to follow the first-order MIT to its 
second order critical end point in the absence of mag- 
netic long-range order. Only recently Mott criticality 
was also investigated for the spin degree of freedom by 
nuclear magnetic resonance (NMR) measurements under 



pressure.— The focus of the latter study was on the critical 
scaling behavior upon varying pressure at fixed temper- 
ature, which is described by a critical exponent S. Here 
we present a theoretical modeling of the spin and charge 
Mott criticality and determine 6. The spin criticality will 
be investigated via the local spin susceptibility, which is 
related directly to the NMR spin-lattice relaxation rate 
1/Ti. For the charge criticality we focus on the local 
density of states at the Fermi energy. 

Experimentally Mott criticality is probed by the scal- 
ing behavior of a selected quantity a (e.g., the conduc- 
tivity) as a function of external parameters such as the 
temperature T or the pressure p near the critical end 
point. Specifically one observes scaling with respect to 
the reduced parameters tred = (T — Tc)/Tc and Pred = 
{p — Pc)/pc, where the index c denotes the values at the 
critical end point. Criticality is then classified by the set 
of exponents /3, 7, and 6 via 
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where ctc = cr(ti.cd = O,prod = 0). The critical exponents 
obey the scaling law 7 = J(/3 — 1)4 

The Mott transition has been proposed to be in the 
Ising universality class^"— based on the assumption that 
the double occupancy may serve as a fingerprint observ- 
able for the MIT, which plays a similar role as the or- 
der parameter in a thermodynamic phase transition to 
a broken symmetry state. The universality class of a 
second-order phase transition is determined by the sym- 
metry of the order parameter and the spatial dimension; 
the scalar character of the double occupancy would imply 
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the Ising universahty classi^ Conductivity measurements 
for (Vi_2;Cra;)203 indeed confirm critical behavior com- 
patible with 3D Ising universality (/? « 0.33, 7 « 1.2, 5 
w 4.8) but the situation for the two-dimensional k-C1 
has remained controversial. 

The critical exponents of the 2D Ising model are /3 = 
1/8, 7 ~ 7/4, and S = 15J^ Conductivity measurements 
under pressure performed on k-C1 challenge the predic- 
tion of Ising universality, since the observed exponents 
are /3 « 1, 7 w 1, and 5 w 2. Imada et al. argued that 
the observed deviation from Ising universality is a mani- 
festation of unconventional quantum criticality specific to 
a two-dimensional systemJ^ii^ A different scenario was 
proposed by Papanikolaou et al. who claimed that the 
conductivity can have a different critical behavior than 
the order parameter of the transitions^ 

In NMR experiments under pressure on k-CI Ka- 
gawa et al. observed that the critical enhancement of 
the conductivity upon passing through the critical end 
point is accompanied by a critical suppression of spin 
fluctuationSfi the latter was inferred from a decrease of 
the nuclear spin- lattice relaxation rate T^^. Identical 
critical exponents d were determined for the conductiv- 
ity and the spin relaxation \/{TiT) within experimental 
accuracy. 

Here we aim at a microscopic description of Mott 
criticality in organic conductors by studying the two- 
dimensional one-band Hubbard model on an anisotropic 
triangular lattice with the Hamiltonian 



(2) 



where ^ {^k a) creates (annihilates) an electron in a 
Bloch state with lattice momentum k. is the local 
density operator for site i and spin a ~ "[.i^.U > Q is the 
local Coulomb repulsion strength, and /x is the chemical 
potential. The electronic dispersion is given by 



Efc — — 2t (cos kx + cos ky) — 2tdiag c,oii{kx + ky 



(3) 



Following Ref. [la we choose for k-C1 a diagonal hopping 
^diag = 0.44i and fix the filling at n = 1 in a grand- 
canonical calculation. 

We obtain an approximate solution of the Hubbard 
model by using cluster dynamical mean-field theory 
(CDMFT) on a 2 X 2 plaquette. The CDMFT self- 



consistency equationaiiii^ are 



G{iuJn) = ^ ((iwn + fJ.)l- Hiiujn) ~ t{k) 

k 



^o^(iwn) = G ^(iw„) - i:(iw„ 



(4) 
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For the = 2 X 2 plaquette CDMFT, the hop_ping 
matrix t{k) is defined via its matrix elements tij{k) = 



,i(fc+fc)-(X,-X 



e 



fe-Hfc' 



r, where Xi and Xj are the 



N.-'Ek 

position vectors of cluster sites i and j, k is in the re- 
duced Brillouin zone, and the cluster momenta take the 



values k = (0,0), (7r,0), (0,7r), and (tt, tt). All quantities, 
i.e. t, the coarse-grained cluster Green function G, the 
Weiss field Qo, and the cluster self-energy S are Nc x Nc 
matrices, and 1 is the unit matrix. In the following we 
consider only paramagnetic solutions and the spin index 
is therefore suppressed. 

The self-consistency cycle is closed by solving the im- 
purity problem, i.e. by calculating a new cluster Green 
function matrix Gij{T) — ~{7'tC.i{t)c^j{0))s^ii for a given 
self-energy and the corresponding Weiss field. Here S'cff 
denotes the effective action of the auxiliary Anderson 
impurity model, which is solved by numerically exact 
continuous-time quantum Monte Carlo (QMC) simula- 
tions based on the expansion of SeS in the impurity-bath 
hybridization^^"— In contrast to single-site DMFT;^"— 
CDMFT takes short-ranged nonlocal correlations within 
the cluster into account. These nonlocal correlations are 
particularly important for two-dimensional systems 

We employ the following strategy for obtaining infor- 
mation on the critical behavior at the continuous Mott 
transition: First, we calculate the double occupancy D 
— Sil^^i ^ function of U for a fixed ratio 

t/T and search for hysteresis, i.e. whether there is a finite 
U region in which both a metallic and a Mott-insulating 
solution of the self-consistent CDMFT equations exist de- 
pending on the initial guess for the self-energy. If hystere- 
sis occurs, the temperature T is increased, otherwise T is 
decreased. This procedure is repeated until the boundary 
between hysteretic and non-hysteretic behavior is deter- 
mined accurately. The critical inverse temperature for 
the continuous Mott transition is denoted as {t/T)c\ the 
critical end point is determined by the two parameters 
{U/T)c and {t/T)c. Spin and charge criticality are sub- 
sequently measured and quantified by 



Gioc(i/(2r)), 

Im Xioc(w) 



lim 



(6) 
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which are both functions of the reduced variables tjod and 

Prcd- 

The local Green function Gioc(r) = N'^ E!I'i(ci(T-)4> 
measured at imaginary time r 1/{2T) approximates 
TA{u! = 0) and thus gives an estimate for the local den- 
sity of states at the Fermi energy without necessitating 
an analytical continuation procedure for the imaginary- 
time datai^i^ Therefore (Jch serves as one possible mea- 
sure for the criticality of the charge degree of freedom. 
The spin excitation spectrum is refiected in the local dy- 
namical spin susceptibility xioc {^) , which is calculated by 
a QMC measurement of the imaginary time correlation 
function xioc(t) = N^^ Y.f=i{Si,z{T)Si^z{0)) and the an- 
alytic continuation of its Matsubara transform to real fre- 
quencies. Here we use the maximum entropy method^ 
for the bosonic kernel according to 



X1oc(t) 



duj e 



-lmxioc(w)- (8) 



In the Mott insulator Imxioc('^)/w is sharply peaked 
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T/U 



(T/U),^l/52 



critical 

path through the continuous 
end point Mott transition at fixed 
temperature 




rf/r;,=7.68 



FIG. 1. Schematic path through the continuous Mott tran- 
sition at fixed temperature. We assume that the pressure- 
controlled transition in Ref. corresponds approximately to 
the bandwidth-controlled transition in the Hubbard model. 
The hopping amplitude t (bandwidth W oc t) is varied (blue 
arrow) while the local Coulomb repulsion U, the relative 
anisotropic diagonal hopping tdiag/t and the temperature T 
are kept fixed in our calculations. The indicated values of 
{T/U)c and {t/T)c are the critical values found for tdiag/t ~ 
0.44. 

at a; = 0. This behavior prohibits a reliable determina- 
tion of (7sp in the insulator from our numerical data. We 
therefore restrict the analysis of spin criticality to the 
metallic side of the transition. 

In order to model the bandwidth-controlled Mott tran- 
sition by tuning the pressure, some further assumptions 
are necessary: We assume that varying t/U amounts to 
varying pressure, and that T/U is kept fixed at constant 
temperature. In essence this implies that the value of U 
is fixed independent of external conditions in the exper- 
iment. This is motivated by the fact that the Hubbard 
interaction is strictly local. Moreover we make the ap- 
proximation that tdiag/t = 0.44 remains fixed at the value 
taken from a fit to the band structure^^ even when pres- 
sure is applied. Fig. [T] summarizes our strategy for mod- 
eling the continuous Mott transition across the critical 
end point. 

The basis for the discussion of the critical behavior are 
the data displayed in Fig. [5] Both charge and spin de- 
grees of freedom show critical behavior with an infinite 
slope at the continuous transition, which is identified at 
{t/T)c = 7.68, independently for charge and spin. The 
theoretical results resemble the experimental data for the 
criticality of the conductivity and the NMR spin-lattice 
relaxation rate 1/Ti in Ref. qualitatively. The increase 
of Gioc(l/(2r)) indicates that the low-energy spectral 
weight increases upon passing from the insulator to the 
metal. This behavior reflects the closing and the filling of 
the charge gap with low-energy states at finite tempera- 
tures. Similarly, the measured conductivity will increase 
in the metal. In contrast, the spin susceptibility behaves 
oppositely; it is suppressed with increasing "pressure" on 
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FIG. 2. Critical behavior at the continuous bandwidth- 
controlled Mott transition (vertical dashed line, {T/t)c = 
7.68). Upper panel: Evolution of [och — Oc.ch)/uc,ch and 
(o"sp — o-c,3p)/o'c,3p upon increasing t/T. The solid curves are 
fits to the critical scaling behavior for the exponent S accord- 
ing to Eq. ([l}. Lower panel: Double-logarithmic plot of the 
same data measured from the critical end point of the Mott 
transition. Solid lines show the scaling fits (the same fits as 
in the upper panel). Error bars are estimated from the linear 
regression fit to the data in the double-logarithmic scale. 

the metallic side of the continuous MIT and enhanced in 
the insulator, as measured by l/(riT). 

This qualitative behavior of the spin susceptibility 
finds a natural interpretation in terms of the probabil- 
ities of relevant plaquette cigenstatcs in the ensemblci^ 
The enhancement of spin fluctuations in the insulator is 
thereby traced to the predominant occupation of the pla- 
quette with a four-electron singlet stated with zero total 
momentum and zero total spin S = 0. The second-most 
probable states are the three triplet states with spin S 
= 1. Since the singlet state has a high occupation prob- 
ability in the insulator, spin flip {A S = 1) excitations 
to the triplet states are more likely and lead to the large 
susceptibility in the insulator. 
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TABLE I. Summary of critical exponent S for various models 
and experiments. 
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The lower panel of Fig. [2] shows the critical scaling 
behavior in a double-logarithmic plot of a^h and a^p rel- 
ative to their values at the critical point as a function 
of \t/T - {t/T)^\/{t/T)c. From linear fits to the double- 
logarithmic plot we extract the critical exponents accord- 
ing to Eq. ([T]). Error bars are estimated from the linear 
regressioir fit to the data in the double-logarithmic scale. 
For the charge criticality the exponent 1/5 = 0.58 ± 0.06 
is obtained on the insulating and 1/5 = 0.61 ± 0.05 on 
the metallic side of the transition. The critical behav- 
ior with an infinite slope of the density of states at the 
Fermi energy is apparent and resembles the measured dc 
conductivity!^ The virtue of a direct quantitative com- 
parison is, however, uncertain because the conductivity 
and the density of states at the Fermi energy may follow 
different scaling lawsj^ For the spin criticality, instead, a 
comparison is meaningful and the quantitative agreement 



with experiment is remarkable, with an exponent 0.48 ± 
0.03 consistent within error bars with the experimentally 
determined 1/5 « 0.5.^ 

Table |T] summarizes the values of the critical exponent 
(5 > 1 for various theoretical models and for several ex- 
periments probing the continuous Mott transition. Note 
that (5 — > 1 corresponds to a smooth transition with a van- 
ishing discontinuity in the first derivative of CTch, while 5 
— > oo describes a discontinuous transition with a discon- 
tinuity in fjch itself. 

Our results for plaquette CDMFT are consistent with 
the experimental data for re-Cl both for the charge and 
the spin degree of freedom. In particular, the critical 
exponent 5 is much smaller than expected for the Ising 
universality class in two dimensions [5 ~ 15) and also 
smaller than in single-site DMFT (5 ^ 3). The dif- 
ference between single-site DMFT and four-site cluster 
DMFT is indeed striking. The latter includes nonlocal 
correlations and thereby allows for a coarse momentum- 
space differentiation. Hence, electrons in parts of the 
Brillouin zone become localized already in the metallic 
phase. In contrast to single-site DMFT, where the Mott 
transition occurs simultaneously by a loss of quasipar- 
ticle integrity along the whole Fermi surface, the Fermi 
surface itself may disintegrate and ultimately vanish at 
the Mott transition within cluster DMFT,^ Therefore, 
only a smaller fraction of electrons indeed exhibits crit- 
ical behavior, which is why the transition is smoother 
(translating into a smaller value of 5) than in single-site 
DMFT. 
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